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Abstract 

Purpose Integrating soil quality impacts in life cycle assess¬ 
ment (LCA) requires a global approach to assess impacts on 
soil quality that can be adapted to individual soil and climate 
contexts. We have developed a framework for quantifying 
indicators of impact on soil quality, valid for all soil and 
climate conditions, and considering both on-site and off-site 
agricultural soils. Herein, we present one of the framework’s 
impact indicators, which has not yet been quantified in detail 
in LCA studies: soil compaction. 

Material and methods The method includes guidelines and 
tools for estimating midpoint compaction impacts in topsoil 
and subsoil as a loss of soil pore volume (in cubic metre per 
functional unit). The life cycle inventory (LCI) and life 
cycle impact assessment are based on simulation modelling, 
using models simple enough for use by non-experts, general 
enough to be parameterised with available data at a global 
scale and already validated. Data must be as site specific and 
accurate as possible, but if measured data are missing, the 
method has a standardised framework of rules and recom¬ 
mendations for estimating or finding them. The main model 
used, COMPSOIL, predicts compaction due to agricultural 
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traffic. Results are illustrated using a case study involving 
several crops in different soil and climate conditions: a 
representative pig feed produced in Brittany, France. 
Results and discussion Predicted compaction impacts re¬ 
sult from the combination of site-specific soil, climate 
and management characteristics. The data necessary to 
the LCI are readily available from free soil and climate 
databases and research online. Results are consistent with 
compaction observed in the field. Within a soil type, 
predictions are most sensitive to initial bulk density and 
soil water content. 

Conclusions The method lays the foundation for possible 
improvement by refining estimates of initial soil conditions 
or adding models that are simple and robust enough to in¬ 
crease the method’s capacity and accuracy. The soil compac¬ 
tion indicator can be used in LCAs of bio-based materials and 
of waste management stages that consider composting. The 
framework includes other operational indicators (i.e. water 
erosion, soil organic matter change) to assess impact on soil 
quality. They complement other impact categories, providing 
increased ability to identify “impact swapping”. 

Keywords Compaction -Indicators -Life cycle assessment • 
Soil quality 

1 Introduction 

Soils are an essential component of the global ecosystem, as 
important as the hydrosphere and atmosphere for ecosystem 
functioning, that are subject to a series of threats in the form 
of degradative processes (EU 2006). The sustainable devel¬ 
opment of human activities, in particular agriculture and the 
production of bio-based materials (e.g. food, biofuels, wood, 
textile, paper), requires the maintenance of soil quality. Soil 
quality can be defined by a soil’s capacity to function (Karlen 
et al. 1997) and/or its fitness for a given use (Larson and Pierce 


4^ Springer 






Int J Life Cycle Assess (2013) 18:1316-1324 


1317 


1994; Letey et al. 2003), both of which are influenced by 
activities that directly or indirectly impact soil properties. 

Soil compaction, a major process of soil degradation 
around the world (Oldeman et al. 1991), is the decrease in 
soil volume, in particular the air-filled fraction (Hillel 1998). 
It alters the spatial arrangement, size and shape of clods and 
aggregates and, consequently, the pore spaces both inside 
and between them (Defossez and Richard 2002). It can be 
expressed by bulk density, pore volume, porosity, void ratio 
or soil strength. Compaction is one of the eight threats to 
soil identified in the Thematic Strategy for Soil Protection in 
Europe (Huber et al. 2007) but also concerns soils on other 
continents (reviewed by Hamza and Anderson 2005). 
Compaction is found in all systems of land use in which 
land is crossed by machines, people, livestock or wildlife 
and is a major concern in agriculture (Batey 2009). The 
degree of compaction depends not only on the amount and 
surface area of pressure applied, but on the water content 
and weight-bearing capacity of the soil at the time. 
Compaction alters soil structure, limits water and air infil¬ 
tration and impedes root penetration. Ultimately, compac¬ 
tion can influence soil chemical and biological processes 
(e.g. N and C cycles, greenhouse gas emissions), plant 
growth and yield and soil biodiversity. Topsoil compaction 
in hilly landscapes increases runoff and may increase ero¬ 
sion, with subsequent on- and off-farm environmental im¬ 
pacts. It is a cumulative process; additional machine passes, 
for example, may increase compaction when soil water 
content is high (Bakker and Davis 1995). Since it is rela¬ 
tively reversible with ploughing, topsoil compaction is re¬ 
lated more to shorter term economic and environmental 
impacts, whereas subsoil compaction, difficult to reverse 
even with subsoiling, concerns longer term sustainability. 
In the context of this work, we accept the following defini¬ 
tion for soil: “the naturally occurring, unconsolidated mineral 
or organic material at least 10 cm thick that occurs at the 
earth’s surface and is capable of supporting plant growth” 
(SCWG 1998). This definition excludes artificially displaced 
materials (e.g. mine spoils) but includes those subject to soil¬ 
forming processes. Nonsoil is the superficial materials that do 
not meet the preceding definition of soil. 

Until recently, impacts of production systems on soil 
quality have not been explicitly included in life cycle as¬ 
sessment (LCA) (Brandao et al. 2011) despite recommen¬ 
dations to do so (Mila i Canals et al. 2007a; 2007b). 
Although a scarce resource, soil has been treated as an 
“ancillary” item when estimating land use impacts in LCA, 
contributing to the maintenance of production processes but 
not part of them (Fava et al. 1990). When impacts on soil 
become the subject of study, the soil must be included 
within LCA system boundaries, which means considering 
the state of the soil, not only its inputs and outputs, 
when estimating impacts. Consideration of system state 


already exists in LCA for indicators of impacts on water 
resources (Pfister et al. 2011) and biodiversity (Curran et 
al. 2011). 

In a previous paper (Garrigues et al. 2012), we addressed 
the challenges of integrating an indicator of impacts on soil 
quality into the methodological structure of LCA. Impacts 
on soil quality in LCA are more easily expressed by pro¬ 
cesses than by physical, chemical or biological properties 
and functions because processes are easier to relate to func¬ 
tional units. Erosion and soil organic matter (SOM) impacts 
already exist in LCA approaches (Mila i Canals et al. 2007a; 
Nunez et al. 2010), but compaction impacts have received 
relatively little attention. Cowell and Clift (2000) proposed a 
“soil compaction indicator” for agricultural subsoil that 
equalled the product of the mass of field machinery and 
the time they spend in the field (i.e. potential impact in 
tonne-hours). Although the indicator did not consider the 
compactability of soil, the authors mentioned that it could be 
included (analogous to a fate factor for toxicity impacts). 
More recently, Oberholzer et al. (2012) proposed the Swiss 
Agricultural Life Cycle Assessment - Soil Quality method 
(SALCA-SQ) for including on-farm soil quality in LCAs. 
One of SALCA-SQ’s nine qualitative indicators assesses the 
influence of a given cropping system on soil macropore 
volume. It is a function, among other factors, of subsoil 
compaction risk, itself a function of the greatest single¬ 
axle pressure exerted by each machine (at a depth of 35 
cm), its area of contact, soil firmness and soil water content 
on operation dates (Oberholzer et al. 2006). Compaction 
risk scores (from 0 to -2) from all machines and plots are 
added together to estimate farm-scale compaction risk. Its 
predictions of soil quality impacts of experimental organic 
and conventional crop rotations appear largely consistent 
with field observations (Oberholzer et al. 2012). SALCA- 
SQ, however, is designed to consider soil quality impacts of 
entire cropping systems, not those of the products they 
produce or of products containing several crop-based in¬ 
gredients (potentially from different farms around the 
world). For it to do so would require a significant amount 
of work, since SALCA-SQ is highly oriented and calibrated 
for Swiss conditions (e.g. in its lookup tables). In addition, 
its qualitative scale of impact (—, -, 0, +, ++) has the 
potential to render it relatively insensitive to differing im¬ 
pacts of agricultural practices. 

To address some of these issues, we developed a method 
to quantify soil compaction impacts of agricultural products 
in an LCA context that aims to be valid for all soils and 
climates that considers both on-site and off-site (upstream) 
agricultural soils. It brings together existing water balance 
and soil compaction models in a way that can be reasonably 
applied in LCA context. It takes into account the types and 
dates of field operations, machine mass and tyre width, the 
surface area affected and soil characteristics. Data 
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requirements are modest, allowing the method to be applied 
to products with multiple agricultural ingredients. We ap¬ 
plied the method to a case study of a feed produced for pigs 
raised in Brittany, France. 

2 Methods 

2.1 General framework of the method 

Soil quality impacts are quantified with midpoint impact 
indicators describing processes that can degrade or improve 
the soil. Pathways were selected to link elementary flows of 
the inventory (LCI) to the midpoint indicators, which result 
from the combination of soil, climate and management 
characteristics (Fig. 1). The LCI and life cycle impact 
assessment (LCIA) are based on simulation modelling, 
using models simple enough for use by non-experts, gen¬ 
eral enough to be parameterised with available data for any 
location on Earth and already tested and validated for 
certain scenarios and soils. Input data must be as site spe¬ 
cific and accurate as possible, but if measured data are 
missing, the method has a standardised framework of rules 
and recommendations for estimating or finding them. Most 
of the input data necessary for establishing the LCI are 
common to the first three midpoint indicators (i.e. water 
erosion, SOM change and compaction) (Corson 2012). For 
each indicator, total impact is estimated by summing the 


impacts from individual upstream agricultural sites together. 
Thus, the method currently has no regionalised characterisa¬ 
tion factors for the LCIA, assuming that a given degree of 
erosion, SOM or compaction has equal impact regardless of 
location. For each crop, the temporal boundary includes the 
inter-crop period (if any) that occurs just before the crop. For 
perennial crops, impacts of planting are divided over the 
productive lifetime of the crop to obtain a dimension of 1 year. 
The system boundary for crop products used as raw ingredi¬ 
ents is set at the farm gate, while that for products made from 
these ingredients is set at the factory gate. 

2.2 Soil compaction indicator 

The process to estimate compaction impact for an agricul¬ 
tural production system requires several steps (Fig. 2). First, 
data on machine operations (e.g. operation types and dates, 
machine mass and tyre dimensions) for 1 ha and 1 year of 
crop production are acquired (Electronic Supplementary 
Material, Table 1). Next, soil water content on the dates of 
field operations is estimated. Soil water models require 
knowledge of soil properties that vary spatially and tempo¬ 
rally and are difficult to measure. Thus, the method uses a 
simple water balance model with two reservoirs (BILHY for 
BILan HYdrique, Electronic Supplementary Material, 
Table 2) to calculate the total soil water content to a depth 
of 50 cm (Jacquart and Choisnel 1995). The predicted soil 
water content is assumed to be constant from 0 to 50 cm, as 


Fig. 1 Conceptual steps for 
assessing impacts on soil 
quality (< outlined) (adapted from 
Garrigues et al. 2012). Not all 
potential links between 
midpoint indicators are shown 
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Fig. 2 Operational procedure 
for estimating compaction 
impact (loss of pore volume). 
Asterisk : For fertiliser and 
pesticide applications, each 
subsequent pass begins with the 
final bulk density profile 
predicted after the previous pass 
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the consideration of different horizons (i.e. topsoil vs. sub¬ 
soil) did not increase the accuracy of predictions of soil water 
content in a study by A1 Majou et al. (2008). Input data 
include daily temperature, precipitation and solar radiation, 
which can be acquired from international climate databases 
(NASA 2012). For greater precision, the time step of BILHY 
was decreased from 10 days to 1 day. BILHY requires 
estimates of potential evapotranspiration, which we estimat¬ 
ed with the TURC method (Federer 1996). Then, initial soil 
bulk density is estimated as a function of soil texture. Initial 
dry bulk density comes from the SOTWIS database (ISRIC 
2012), from which soil texture is divided into five classes 
(coarse to very fine, according to the FAO texture triangle), 
each associated with an initial bulk density (Table 1). 

Next, the COMPSOIL model (O’Sullivan et al. 1999; 
Defossez et al. 2003) is used to predict bulk density after each 
machine operation using readily available mass and tyre data. 
The depth to which compressive forces are transmitted de¬ 
pends on soil moisture, bulk density and machine character¬ 
istics. For simplicity, the method assumes uniform initial bulk 
density and water content profiles. Simulations begin with the 
same initial bulk density for all field operations, which implies 
that machines do not follow the same tracks each time. An 
exception is made, however, for fertiliser and pesticide appli¬ 
cations, which tend to follow the same tracks; in this case, 
each subsequent pass begins with the final bulk density profile 
predicted after the previous pass. The passage of ploughs is 


Table 1 Mean dry bulk density (p s ) of texture classes of the FAO 
texture triangle calculated from the SOTWIS database (ISRIC 2012). 
Mean p s is used as initial bulk density in the COMPSOIL model 


Texture 

Coarse 

Medium 

Medium fine 

Fine 

Very fine 

Mean p s 

1.41 

1.35 

1.29 

1.30 

1.25 

n 

474 

1,625 

182 

3,011 

1,519 


assumed to result in no increase in soil compaction from 0 to 
30 cm; only predicted compaction below 30 cm is considered. 
COMPSOIL predicts final bulk density every 5 cm from 0 to 
50 cm deep. Mean topsoil (0-30 cm) and subsoil (30-50 cm) 
final bulk densities are calculated separately. 

The increase in mean dry bulk density is converted into a 
loss of pore volume (A/in cubic metre): 

Af = (/f -/) x As x D, (1) 


where ff and f are the final and initial porosity, respectively 
(in percent), A aff the affected area (in square metre) and D 
soil depth (in metre) 


/ 


Ps ~ ffb _ J _Pb 

Ps Ps, 



Pb = 2.65 x (1 —/), 



where is the dry bulk density, p s particle density=2.65 t/m 

and/the porosity (in percent). 


A 


aff — Av max 


X 




where T w max is the width of the widest tyre of the machine 
operation system (in metre) and W w the fieldwork width 
(in metre), which is needed to calculate the number of 100-m 
passes in a square 1 -ha field. 

For each operation, topsoil and subsoil losses of pore 
volume are added together to equal the total loss of pore 
volume for the operation. Finally, the losses of pore volume 
of all operations are added together to obtain the total loss per 
hectare of crop, which is the value used for LCI and LCLA. 
Loss of pore volume is estimated for crop-based products by 
dividing the impact per hectare by the dry matter yield per 
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hectare. Loss of pore volume for co-products made from raw 
crop products (e.g. soya oil and meal) is calculated as usual, 
allocating impacts to co-products when necessary. 

2.3 Case study 

The case study was selected to illustrate impacts of a compos¬ 
ite product formed from crop-based ingredients produced in 
and with widely differing soils, climates and crop manage¬ 
ment practices. It focused on the global soil quality impacts of 
producing feed for pigs raised in Brittany, France, with in¬ 
gredients (maize, wheat, triticale, barley, pea, rapeseed meal, 
soya meal and oil and cane molasses) coming from Brittany, 
Brazil and Pakistan (Table 2). Since the method requires data 
as precise as possible, we obtained regional data for agri¬ 
cultural practices in France (Brittany) and Brazil (Santa 
Catarina) using an INRA LCA database (UMR SAS, 
Rennes, France) and commercial and online sources (details 
in Electronic Supplementary Material, Table 3). We assumed 
typical management practices and crop yields equal to 5-year 
means. For sugarcane, harvested annually but planted once 
every 6 years, the impacts of the tillage and sowing oper¬ 
ations needed to establish it were divided by six. Soil data 
came from the Harmonized World Soil Database (FAO et al. 
2009). Climate data came from Meteo France for France and 
EPAGRI (Empresa de Pesquisa Agropecuaria e Extensao 
Rural de Santa Catarina) for Brazil (Vamilson Prudencio da 
Silva, personal communication). Pakistani climate data were 
acquired from international climate databases (NASA 2012) 
for coordinates corresponding to the centre of the Punjab 
region (30°N, 70°E). For Pakistan, soil and climate character¬ 
istics were chosen from the centre of Punjab, the country’s 
main production zone of sugarcane. Economic allocation 
(with prices averaged over 2004-2007) was used to allocate 
impacts of canola, soya and sugarcane to canola meal, soya 
meal and oil and molasses, respectively. Impacts per tonne of 
ingredient in pig feed were added together to calculate total 
impact per tonne of pig feed produced in Brittany. 

3 Results 

The total area affected per hectare by compaction was similar 
for crops grown in Brittany (19,620-20,690 m /ha) due to 
similar tillage practices (see Table 2). In the topsoil, maize 
production caused the largest loss of pore volume per hectare 
(185.2 m / ha). For crops in loamy soils, topsoil losses of pore 
volume per hectare tended to be three times those in the 
subsoil, except for canola (topsoiFsubsoil ratio of 3.4 instead 
of 3.0). For soya, the only crop in a clayey soil, the ratio was 
lower (2.4). The range of pore volume losses per hectare 
tended to be small for the crops grown in Brittany (171- 
185 m /ha), and the pore volume losses per tonne of crop 


were inversely proportional to yields. For pig feed ingredients, 
the relative impact of each ingredient tended to be proportion¬ 
al to its relative mass in the feed. 

Predicted impacts in Brittany, France, were consistent with 
the machinery used there for each crop type. For example, 
maize operations need heavier tractors with wider wheels than 
those for winter cereals, but they pass less often over the field. 
In Santa Catarina State, Brazil, reduced tillage is used, de¬ 
creasing the number of machine operations and therefore, 
potentially, the impact. In Brittany, field operation dates for 
maize have a high probability of occurring during periods 
when its loamy soils have relatively high water content. In 
contrast, the clayey soils underneath Brazilian soya compact 
less, even though machines used for soya are heavier than 
those used for maize. In Pakistan, planting sugarcane only 
once every 6 years decreases the compaction impact of plant¬ 
ing over the life cycle of the crop. Of the six terrestrial 
operations considered for sugarcane cultivation (Electronic 
Supplementary Material, Table 3), the impacts of the first five 
are divided by six (crop establishment), three are not consid¬ 
ered because they are made by plane. For rapeseed, operation 
dates correspond to periods of drier soil in Brittany, decreasing 
stress propagation, and thus compaction, in the soil. 

4 Sensitivity analysis 

A “one-factor-at-a-time” sensitivity analysis was performed 
to estimate the relative sensitivity of predictions of pore 
volume loss per hectare due to production of a single crop 
(maize) to individual variations in soil type (loam vs. clay), 
initial bulk density and gravimetric water content. 
Sensitivity was measured with a sensitivity index, calculated 
by dividing the normalised range of output values by the 
normalised range of the associated input values. The sensi¬ 
tivity analysis showed that loss of pore volume was more 
sensitive to initial bulk density and water content in loamy 
soils than clayey soils. Overall, sensitivity to initial bulk 
density increases as soil water content increases, and sensi¬ 
tivity to soil water content increases as initial bulk density 
decreases (Fig. 3). Future sensitivity analyses can be done 
for other factors, such as machine mass and width, tillage 
practices for different type of location and humidity. 

5 Discussion 

The compaction indicator is currently restricted to life cycle 
stages of agricultural production that involves machine field 
operations. The trends in comparisons of potential compac¬ 
tion estimates are consistent with general experimental ob¬ 
servations: compaction due to traffic increases with soil 
moisture up to a specific water content and depends on the 
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Fig. 3 Sensitivity of loss of 
pore volume (in cubic metre) in 
soil due to machine traffic to 
produce grain maize in Brittany, 
France, to initial soil bulk 
density (in tonnes per cubic 
metre) and gravimetric water 
content (in percent) for a loam 
or b clay soil. Note that 
maximum soil water content 
changes with soil type 
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mass and tyre width of the machinery used (Richard et al. 
1999), soil type, tillage management and climate (Boizard et 
al. 2001). Although the COMPSOIL model has been vali¬ 
dated only for moderate loads on sandy loam and clay loam 
soils (O’Sullivan et al. 1999), its function within the method 
is to identify any relative differences in compaction impact 
between scenarios, not to make highly accurate bulk density 
predictions. The sample case study, quantifying pore vol¬ 
ume loss due to production of different crops in different 
countries, shows the feasibility of the method at different 
levels of data availability (high for France, medium for 
Brazil and low for Pakistan). 

Many of the necessary input data (e.g. agricultural prac¬ 
tices, field operations, operation dates) are standard input 
data for crop LCAs. Data about agricultural machines (e.g. 
mass, tyre width, operating width) may be more difficult to 
find, but many can be found online. The data collected for 
this research form one such database. We believe that the 
method reaches a balance between being overly complex 


and overly simple. Simplifications include assuming uni¬ 
form water content and bulk density profiles and, for most 
operations in a sequence, the same initial bulk density. The 
uniform bulk density profile may not represent physical 
reality, but this simplification seems acceptable given un¬ 
certainty and lack of knowledge about changes in bulk 
density with depth and the fact that the indicator aggregates 
bulk density changes throughout the profile. The use of 
COMPSOIL model depends on the soil type. The model 
requires a stress/strain relationship and mechanical soil pa¬ 
rameters and their variations as a function of a variety of 
physical properties. The constants for estimating soil param¬ 
eters of model equations are given by O’Sullivan et al. 
(1999) for sandy loam and clay loam and by Defossez et 
al. (2003) for loess and calcareous soil. Saffih-Hdadi et al. 
(2009) gives pedotransfer functions of representative soils 
covering the soil texture triangle. 

In addition, BILHY, developed to identify soil water defi¬ 
cits, cannot simulate water contents greater than field capacity, 


4^ Springer 







Int J Life Cycle Assess (2013) 18:1316-1324 


1323 


such as saturated soil. Replacing the two-compartment model 
of soil water content with a two-dimensional soil model could 
predict the soil water profile more precisely (i.e. every 5 cm), 
but the added complexity may not result in greater accuracy. 
Studies indicate that an increase in model realism may not 
only fail to increase model accuracy, but may decrease it (and 
model generality, as well) (Weisberg 2006). Water content in 
the first few centimetres of soil may greatly influence rutting 
and compaction impacts because that is where wheels contact 
the soil, but the extra parameters required to simulate two- 
dimensional water flow, and the uncertainty around them, do 
not currently warrant burdening the method with further input 
data requirements. Nonetheless, the compaction indicator cur¬ 
rently considers only a cumulative and negative impact. The 
indicator could be expanded to consider a positive impact of 
ploughing on soil pore volume, balanced with a negative 
impact on long-term soil structure. The compaction indicator 
could also be expanded by considering compaction due to 
livestock trampling and the impact of climate on compaction 
resilience and on bare soil. The long-term cropping impact 
could be considered, for example, by setting initial soil con¬ 
ditions as a function of previous tillage practices or consider¬ 
ing the impact of alternating root systems of different depths. 
Each of these improvements would require additional re¬ 
search. COMPSOIL does not explicitly consider the influence 
of SOM content on compaction, but the framework (Garrigues 
et al. 2012) already includes indicators of SOM change and 
erosion to assess impacts on soil quality (Corson 2012). All 
three indicators can be used in LCAs of bio-based materials, 
for example to compare cropping or waste management tech¬ 
niques, such as composting. Although the impact of non¬ 
cultivation processes on soil quality could be included, many 
of them, such as soil sealing with concrete, transform soil into 
nonsoil, which has zero soil quality. Thus, we believe that 
non-agricultural processes are better included with impacts of 
land use and land use change. The soil quality impact in¬ 
dicators can interact with other impact categories, such as 
climate change, in which SOM changes influence net C emis¬ 
sions into the atmosphere, and eutrophication and biotic pro¬ 
duction, which are influenced by erosion. Also, soil quality 
impact indicators complement other impact categories, pro¬ 
viding increased ability to identify “burden shifting” or trade¬ 
offs between transport distance and soil quality. All the input 
data necessary for establishing the LCI for the three indicators 
(approximately 30 parameters) are presented in the project 
report (Corson 2012). Most input data are common to the 
three indicators, and the inclusion of new soil quality impact 
indicators (such as salinisation) should only slightly increase 
input data requirements. The soil compaction indicator could 
be linked to midpoint/endpoint indicators in the future by 
specifying impact pathways between a change in soil porosity 
and water infiltration, gas exchange with the atmosphere, and 
biological processes. 


6 Conclusions 

Impacts on soil quality should be taken into account into a life 
cycle perspective because of the essential role of soils in 
ecosystem functioning. Integrating soil quality impacts 
throughout the life cycle of an agricultural product requires a 
global approach to assess impacts on soil quality that can be 
adapted to individual soil and climate contexts. We have 
developed a framework for quantifying indicators of impact 
on soil quality, valid for all soil and climate conditions, and 
considering both on-site and off-site agricultural soils. These 
indicators can be used in LCAs of bio-based materials or the 
waste management stage that considers composting. Soil com¬ 
paction has not been quantified in detail in LCA. Total com¬ 
paction impact results from the combination of soil, climate 
and management characteristics and is expressed as a loss of 
pore volume (in cubic metre/FU). Most input data can be 
found in existing national and international databases, except 
for some agricultural practices, which vary widely by region. 

The method aggregates compaction impact of all crop- 
based ingredients in different soil and climate conditions. 
Estimates of compaction impact are consistent with empirical 
observations, increasing with soil water content and depending 
on soil type and machine mass and tyre width. Within a soil 
type, compaction predictions are most sensitive to initial bulk 
density and soil water content. 

The method developed answers needs identified by 
Garrigues et al. (2012) for LCA indicators of impacts on soil 
quality. The framework allows for incremental improvement, 
for example by refining the initial state of the soil as a function 
of previous tillage practices. The method currently considers 
only compaction due to agricultural machines, but other causes 
of compaction could be added (e.g. animal trampling) if the 
addition is simple enough to be used by non-experts, general 
enough to parameterise with available data at a global scale and 
sufficiently validated. 
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